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processes  belong  to  the  class  of  stationary  linear  time  series.  A  linear 

m 

time  series  -»• h(t-s)e(»)J,  is  Gaussian  if  the  distribution  of  the 

independent  innovations  (e(t)}  is  normal.  Assuming  that  Ee(t)  -  0,  some  of 
the  third  order  cumulants  c___(m,n)  ■  Ex(t)x(t+m)x(t+n)  vrt.ll  be  non-zer  if 
the  e(t)  are  not  normal  and  E^^(t)  *  0.  If  the  relationship  between  {x(t)} 
and  {e(t)>  is  non-linear,  then  {x(t)}  is  non-Gauss ian  even  If  the  e(t)  are 
normal.  This  paper  presents  a  simple  estimator  of  the  bispectrum,  the 
Fourier  transform  of  {c__(m,n)}.  This  sample  bispectrum  is  used  to 

XXX 

construct  a  statistic  to  test  whether  the  blspectrum  of  (x(t)}  is  non-zero. 

A  rejection  of  the  null  hypothesis  implies  a  rejection  of  the  hypothesis 
that  (x(t)}  is  Gaussian.  A  related  test  statistic  is  then  presented  for 
testing  the  hypothesis  that  {x(t)}  is  linear.  The  asymptotic  properties  of 
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l.  Introduction 

Tine  series  data  from  a  variety  of  sources  are  often  analyzed  under 
Che  explicit  or  implicit  assumption  that  they  are  generated  by 
autoregressive  (AR)  or  autoregressive  moving  average  (ARMA)  processes. 
These  processes  are  finite  parameter  linear  stationary  stochastic 
processes.  The  tth  element  of  a  causal  linear  processes  {x(t)>  is  of  Che 
form 

m 

x(c)  -  l  h(s)e(t-s) ,  (1.1) 

3-0 

where  the  e(t)  are  independent  Identically  distributed  random  innovations 
with  Ee(t)  -  0.  In  filter  theory  terminology,  the  stationary  process 
(e(t)}  is  Che  input  to  a  time  invariant  linear  filter  whose  impulse 
response  is  (h(t):  t  -  0,  1,...}.  If  Ic»Qh*(t)  <  “,  the  covariance 
function  of  the  stationary  output  process  (x(t)}  is  finite.  If  the  input 
is  Gaussian,  then  the  output  is  Gaussian  and  its  covariance  function 
completely  determines  the  joint  distributions  of  the  process. 

But  suppose  that  the  e(t)  are  not  normal  and  M3  -  Ee^(t)  *  0.  Then 
the  third  order  cumulant  c^xxCm.n)  -  Ex(t)x(t+m)x(t+n)  *  0  for  many 
values  of  m  and  n.  The  same  is  true  if  {x(t)>  is  generated  by  a 
nonlinear  filtering  operation  satisfies  a  Volterra  functional  expansion 
(Brlllinger,  Sec.  2.10,  1975).  Nonlinear  models  are  beginning  to  play  a 


2 


role  in  applied  time  series.  An  overview  of  nonlinear  models  is  given  by 
Priestley  (1980). 

Using  an  estimator  of  the  bispectrum  of  {x(t)>  Subba  Rao  and  Gabr 
(1980)  present  tests  for  whether  the  process  is  Gaussian  and  whether  it 
is  linear.  They  do  not  use  the  asymptotic  variance  of  the  sample 
bispectrum  in  their  multivariate  procedures,  which  is  reasonable  for 
smallish  sample  si2es.  This  paper  presents  a  modification  of  their 
approach  that  malces  heavy  use  of  the  large  sample  properties  of  the 
sample  bispectrum.  Let  us  begin  with  a  brief  review  of  the  bispectrum  of 
a  stationary  zero-mean  process. 

2.  The  Bispectrum 

The  bispectrum  B(wpU2)  gives  a  measure  of  the  multiplicative 
nonlinear  interaction  of  frequency  components  in  {x(t)>  (Hasselman  et 
al. ,  1963).  For  a  real  stationary  time  series  with  Ex(t)  «  0,  the 
bispectrum  is  defined  as  follows: 

m  m 

V“i’“  2)  ■  l  l  cxxx(m,n)exp[-i(utim  +  ui2n)  1 ,  (2.1) 

m«-»  n»-« 

assuming  that  |cxxx(m,n)|  is  sumraable.  Given  the  symmetries  of 
Bx(iitpU2),  its  principle  domain  is  the  triangular  set  fl  »  (0  <  <  ir , 

U2  <  up  +  u>2  <  2tt}  (Van  Ness  1966).  The  area  of  this  triangle  is 
2/3  of  the  area  of  the  triangle 

If  (x(t)}  is  linear,  then 


Bx(up<u2)  -  y3H(u)1)H(<j)2)H*(u1-hu2) 


(2.2) 
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where  H(w)  ■  Jt.(jh(t)exp<-iwC)  is  Che  filter  transfer  function  and  the 
scar  denotes  complex  conjugate  (Brillinger,  1965).  Thus  if  u3*Ee^(t)#0, 
then  Bx(uiiM2^  ^ 

The  finite  Fourier  transform  of  a  sample  (x(0),  x(l) .  x(N-l)}  of 

the  process  can  be  used  to  construct  a  consistent  estimator  of  the 
bispectrum. I  Let  ■  2xn/N  for  n  »  0,  1,...,  N-l.  For  each  pair  of 
integers  j  and  k,  define 

F(j,k)  -  (2.3) 

where 

N-l 

X(u».«)  -  l  x(t)exp(-iui,t). 
t-0 

Since  X(uj+f()  *  X(uj )  and  X(<uN_j)  ■  X*(<dj),  the  principal  domain  of 
F(j  »k)  is  the  triangular  grid  D  -  {0  <  j  <  N/2,  0  <  k  <  j,  2j  +  k  <  N} 
(let  N  be  even) . 

Set  X(0)  -  0,  which  is  equivalent  to  subtracting  the  sample  mean 
from  Che  data.  Thus  F(j,0)  ■  F(0,k)  =  0.  Given  a  summability  condition 
for  the  cumulants  of  (x(t)} ,  it  follows  from  expression  (4.3.15)  in 
Brillinger  (1975)  that 

E(F(J,k)J  -  B*(«j,  Wfc)  +  0(N-1),  (2.4) 

and  the  complex  variance  is 

^Subba  Rao-Gabr  use  the  standard  windowed  sample  covariance  method. 
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EjF(J.k)  -  u»k)|2  -  NSx(uj)S3C(a»k)Sx((Uj+lc)  (2.5) 

(1  +  6(j-k)  +  6(N-2j-k)  +  46(N-3j)6(N-3k)J  +  0(1), 

where  Sx(u)  is  Che  spectrum  of  (x(t)}  and  6(k)  -  0  unless  k  •  0  when 
6(0)  ■  1*  The  complex  covariance  between  F(j  ,k)  and  F(j',k')  is  0(N”2) 
for  j  *  j*  or  k  *  k'  (it  is  0(N“*)  for  the  "exclusive  or"  cases). 

There  are  many  ways  to  average  the  F(j ,k)  to  obtain  a  consistent 
estimate  of  the  bispectrum  on  a  lattice  of  points  in  the  set  D.  For  the 
lattice  L  -  { (2m-l)M/2, (2n-l)M/2:m,n  -  1,2,  ...  and  M  -  N°  for  1/2  <  c  < 
1}  in  D,  a  simple  approach  is  to  average  the  F(j ,k)  in  a  square  of  M2 
points  centered  at  ((2m-l)M/2,  (2n-l)M/2)  if  all  the  ( j ,k)  are  In  the 
domain  (Fig.  1).  The  estimator  is  then 

-  nM-1  nM-1 

BjOu.n)  -  if*  l  l  F(j  ,k)  (2.6) 

j“(m—l)M  k-(n-l)M 

If  a  square  has  points  outside  Che  set  D,  those  points  are  not  Included 
in  the  average. 

If  Bx  (iDj.u^)  is  slowly  varying  over  a  square  of  width  2irM/N  and 
S^w)  is  slowly  varying  over  a  2irM/N  band,  it  follows  from  (2.4)  that 

A 

E[Bx(m,n) ]  -  Bx(2*(2m-1)M/2N,  2*(2n-l)K/2N)  +  0(M/N).  (2.7) 

From  (2.5),  the  complex  variance  is 

A  A 

VarBxU.n)  -  E| Bx(b.,a)  1 2  «  |Bx(».n)|2 

-  NH_4Q{S3t(2x(2m-l)M/2N)  (2.8) 

Sx(2w(2n-l)M/2N)Sx(2x(m+n-l)M/N)  +  0(M/N)] 


where  Q  Is  the  number  of  (j ,k)  in  Che  square  that  are  in  D  but  not  on  the 
boundaries  j  »  k  or  2j+k  ■  N,  plus  twice  the  number  on  these  boundaries. 


If  the  square  is  within  D,  Q  -  M^.  For  any  square,  Q  >  M^/8 
since  the  smallest  set  in  the  domain  (a  triangle)  is  for  the  square 
centered  at  (0,0).  Since  M  -  N°  for  1/2<c<1,  M/N+0  and 

A 

NM  <  NM”^  ■  N*~^c+0  as  N+<».  This  implies  that  B  is  a  consistent 
estimator  of  the  bispectrum  at  (mp^)  as  for  the  sequence  (m(N)  » 

n(N)  »  [u>2n1“CD  .  where  [x]  denotes  the  inceger  part  of  x. 
The  asymptotic  distribution  of  each  estimator  is  complex  normal 
since  che  estimator  is  asymptotically  equal  to  the  one  analyzed  by  Van 
Ness.  Applying  Theorem  4.4.2  in  Brillinger  (1975),  the  estimators  are 
asymptotically  Independent.  Thus  from  (2.8),  the  distribution 
of 

\,n  *  (I,1~4cQ)“1/2[Sx(2iT(2m-l)M/2N)Sx(2ir(2n-l)M/2N) 

Sx(2x(m+n-l)M/N)J'l/2Bx(n,n)  (2.9) 

is  complex  normal  with  unit  variance.  Consequently,  2|Xu,>nl*  is 
approximately  chi-square  with  two  degrees  of  freedom  and  noncentrality 
parameter 

xm,n  "  2(Nl_4cQ)_1Yx(®.a)  >  2N2c-lYx(B,n)  (2.10) 


where 


„  ,  „  2x(2m-i)M  ,  2ir(2n-i)M  .  2*(m+n-l)M  /2ir(2m-l)M  2ir(2n-l)Mx 

Y  (m,n)  -  S  l( - )S  l( - )S"l( - )  B*( - » - ) 

x  2N  x  2N  x  N  ^  2N  2N 
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Moreover,  the  statistic  S  -  2£  |x_  n|2  is  approximately  xiL(X)  where 

(m,n)eL  ’  2P 

X  ■  J  a  an<*  *  the  num*Jer  (®»n)  in  L.  Since  n  »  1, . . . ,  [N/2MJ , 

(a,n)eL  ’ 

P  is  approximately  N^/12M  . 


3.  Testing  for  Gauasianity 

The  statistic  S  is  basically  the  Subba  Rao-Gabr  test  statistic  with 
the  asymptotic  variance-covariance  matrix  instead  of  their  sample 

estimate.  Odder  the  null  hypothesis,  3<a>^,ti>2)  =  0,  and  thus  S  is 
2 

approximately  x  (0)  for  large  N.  This  suggests  the  following  test 
2P 

A  A  A 

statistic:  S  -  2£lttfn|XtB>ll|2  where  ^jfI1  is  given  by  (2.9)  with  S x(u>)  replaced 
by  an  estimate  S^w)  for  each  lattice  frequency  in  the  domain.  If  the 
spectrum  is  estimated  by  averaging  M  adjacent  periodogram  ordinates 

A 

(Fuller,  Sec.  7.2,  1976),  then  Sx(w)  *  Sx(w)[l  +  (M/N)Y]  where  Y  is 
(approximately)  a  standard  normal  variate.  It  then  follows  from  (2.9)  that 

A  A 

S  ■  S  +  Op(M/N) ,  and  thus  the  distribution  of  S  is  also  approximately 
2 

X  (X)  for  large  N.  An  approximate  o-level  test  of  the  null  hypothesis 
2P 

that  B((i»p«*2^  2  0  18  t0  r*J«ct  it  if  S  >  ta  where  o  -  Pr(x2^  * 

If  the  null  hypothesis  is  rejected,  then  the  Gaussian  assumption 
must  be  rejected.  If  not,  then  the  process  may  be  non-Gaussian  but  the 
data  is  consistent  with  a  zero  bispectrum. 

I  will  now  show  that  the  test  is  consistent  as  For  simplicity, 

suppose  that  all  the  squares  are  in  D  so  that  Q  -  M2.  Set 


o  -  (2x)"2//^[Sx(wi)Sx(<02)Sx(u>1-t,«2)  l”1  | ^(<0 1 1 2d<*»  1  *dw2 » 

where  n  is  the  principal  domain.  From  (2.10), 


J 
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P-1  l  Xm>n  -  2N2c-la  +  0(Nc_i) 
(m,n)cL 


(3.2) 


using  Che  incegral  approximation  of  a  sum. 

This  approximation  of  X/P  will  now  be  used  to  obtain  the  large  sample 

property  of  the  power  of  the  chi-square  test.  For  large  N  and  thus  large  P, 
2 

(2P)”1x2P<*)  19  approximately  normal  with  mean  1+(2P)-1\  and  variance 
(2P)  (4P+4X).  Thus  the  large  sample  power  of  the  test  is  a  monotonically 

increasing  function  of  the  "signal-to-noise"  ratio 


P  +  X/2 
(P+-X)1/2 


[P(X/4P)]1/2 
*  (<xN/24)1/2 


(3.3) 


since  P  *  N2^~cVl2.  Thus  the  test  is  consistent.  Moreover,  the  null 
hypothesis  will  be  correctly  rejected  with  probability  near  one  if 
(aN/24)^2  >4.  If  a  ■  1,  for  example,  then  a  sample  size  of  H  >  384  is 
needed  to  achieve  this  somewhat  conservative  bound  for  high  power. 


4.  Testing  for  Linearity 

If  (x(t)}  is  a  linear  process,  then  Sx(u>)  -  |H(w)|2a2  where 

c 

o2  -  Ee2(t).  Thus  from  (2.2), 
e 


Tx(m,n)  S 


-6  2 
ae  u  3 


(4.1) 


for  all  (m,n)cL.  From  (2.10) 


8 


X-P'1  l  (|Vni2  “  l>  (4.2) 

(a,n)cl> 

is  s  consistent  estimator  of  N2c“^yx  under  the  null  hypothesis  that  (x(t)} 
is  linear.  It  then  follows  from  the  results  in  Section  2  that 

•  a 

S  -  N"1/2  l  (IVqI2  -  1  -  »  (4.3) 

(m,n)eL 

is  approximately  normal  N(0>yx/6)  under  the  null  hypothesis.  This  statistic 

A 

with  Sz  instead  of  Sx  can  be  used  to  test  for  linearity.  The  consistency 
of  this  test  follows  from  the  large  sample  analysis  used  to  show  the 
consistency  of  the  S  test. 

5.  Conclusion 

Simple  tests  for  Gaussianlty  and  linearity  of  a  time  series  have  been 
presented.  The  large  sample  variance  and  covariance  of  the  asymptotically 
normal  blspectrua  estimator  are  used  to  simplify  the  Subba  Rao-Gabr  test 
statistics.  The  asymptotic  properties  of  the  statistic  for  testing  for  a 
zero  bispectrum  have  been  presented.  The  power  of  this  test  is  high  when 
aN  is  large. 
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